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Abstract 



Adding a column of numbers produces "carries" along the way. We show that 
random digits produce a pattern of carries with a neat probabilistic description: the 
carries form a one-dependent determinantal point process. This makes it easy to 
answer natural questions: How many carries are typical? Where are they located? 
We show that many further examples, from combinatorics, algebra and group the- 
ory, have essentially the same neat formulae, and that any one-dependent point 
process on the integers is determinantal. The examples give a gentle introduction 
to the emerging fields of one-dependent and determinantal point processes. 

1 Introduction 

Consider the task of adding a single column of digits: 



We have put a dot to the right of a digit whenever the succeeding addition leads to a 
carry. The remainder (mod 10) is written to the right of the dot. Thus 7 -|- 9 = 16 
results in a dot by the 7 and a remainder of 6. At the end, the total is found by adding 
the number of dots (here, five) and appending the final remainder (here, 2) for a total of 
52. The dots are a standard bookkeeping device used (for example) in the Trachtenberg 
system of speed addition [l6i]. 
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How many carries are typical and how are they distributed? Common sense suggests 
that about half the places have carries and that if a carry occurs, it is less likely that 
the next addition gives a carry. Investigating these questions when digits are chosen 
uniformly at random leads to interesting mathematics. 

As additional motivation, look at the remainder column in the example above. Ob- 
serve that there is a carry on the left (and so a dot) if and only if there is a descent on 
the right (a sequence xi,X2, ■ ■ ■ has a descent at z if > Xj+i). Further, if the original 
column of digits is independent and uniformly distributed, so is the remainder column. 
Thus, the results about carries give the distribution of the descent pattern in a random 
sequence of digits. 

[Section "2] derives the distribution theory of carries by elementary arguments. There 
are nice formulae for the chance of any pattern of carries, and the process is one- 
dependent so a variety of probabilistic limit theorems are available. 

[Section "3] introduces determinantal point processes, shows that the carries process 
is determinantal, and illustrates how standard tools for determinantal point processes 
apply to carries. 

[Section "41 reviews the literature on stationary one-dependent processes, shows that 
all of these are determinantal and further that essentially all the neat formulae for carries 
and descents have versions for any stationary one-dependent process. Connections with 
symmetric function theory are developed. A large class of examples arising from the 
work of Polishchuk and Positselski [53] on Koszul algebras is shown to yield many 
natural examples (including the original carries process). 

[Section "Sl gives combinatorial examples: descents in permutations from both uniform 
and non-uniform distributions, the connectivity set of a permutation, and binomial 
posets. 

[Section "HI generalizes from adding numbers to multiplying random elements in a finite 
group. For central extensions (e.g., the quaternions), there is again a carries process 
that is stationary, one-dependent, and determinantal. Thus explicit formulae and limit 
theorems are available. 

[Section "7l contains proofs of some of our more technical results. It shows that any 
one-dependent (possibly non-stationary) point process on the integers is determinantal. 
It also constructs a family of one-dependent determinantal processes (generalizing many 
examples in earlier sections), and computes its correlation kernel. 



2 Probability theory for carries and descents 

Throughout this section we work base b and so with the alphabet B = {0, 1, . . . , 6 — 1}. 
Let Bi,B2,...,B n be a sequence of randomly chosen elements of B (independent and 
identically distributed). There is a descent at i if Bi > -Bj.+i, 1 < i < n — 1. Let Xi be 1 
or as there is a descent at i and D = {i : Xi = 1} he the descent set of Si, ^2, • • • , -B„. 

The following probabilistic facts are elementary. They are stated for carries but 
verification is easier using descents. 

Fact 1. (Single carries) For any i G [n — 1], 
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Thus, when b = 10, the chance of a carry is .45. When b = 2, the chance of a carry is 
.25. For any base, Var(Xj) = | — 

Fact 2. (Runs of carries) For any i and j with l<i<i+j<n, 
P {X, = = • • • = = 1) = f . ^'j / 



J + 1 

Thus, a run of b or more carries is impossible. Further, 

1 / 1 



Gov {X,, X,+i) = E {XiXi+i) - E{X,)E{X,+i) -- ^ 

Fact 3. (Stationary one-dependence) The distribution of {^i}ie[n-i] is stationary: 
for J C [n — 1], i £ [n — 1] with J + i C [n — 1], the distribution of {Xj}j(zj is the same as 
the distribution of {Xj}j,=j^i. Further, the distribution of {^i}ig[n-i] is one-dependent: 
if J C [n—1] has ji,j2 € J =^ \ji — J^l > 1, then {Xj}j^j are jointly independent binary 
random variables with P{Xj = 1) = ^ — The literature on m-dependent random 
variables is extensive. In particular, a classical central limit theorem 40(] shows the 
following: 

Theorem 2.1. For n > 2, the total number of carries T„„i = Xi + • • • + Xn-i has 
mean (n — 1)(^ — ^), variance -[2^(1 — ^) and, normalized by its mean and variance, 
Tn-i has a standard normal limiting distribution for n large. 

Remarks. 

1. An 0(n~^/^) error bound in this central limit theorem can be proved using the 
dependency graph approach to normal approximation by Stein's method 13], 181 ]. 

2. Stationa ry p airwise independent processes can fail to obey the central limit theo- 
rem j3], |43l ]. 

Fact 4. (/c-point correlations) For ^4 C [n — 1], let 

p{A) = P{Xi = 1 for i G ^}. 

For \A\ = k, the p{A) are called A;-point correlations for the point process Xi, . . . , Xn-i- 
They are basic descriptive units for general point processes [17|. For the carries and 
descent process, they are simple to describe. Break ^ C [n — 1] into disjoint, non-empty 
blocks of consecutive integers A = Ai L) A2 U ■ ■ ■ U A^. Thus A = {2, 3, 5, 6, 7, 11} = 
{2, 3} U {5, 6, 7} U {11}. From Facts 2 and 3 above, 

k k 

piA)=fl\( li)A"H a A =(j Ai with \A\=ai. 



i=l 



Fact 5. (Determinant formula) Let ei, £2, . . . , Cn-i be a fixed sequence in {0, 1}. Then 
P {Xi = 61, . . . , = en-i} = ^ • det ("^^^ "ft!: t ^ " ^) • 

Here, if there are exactly k I's in the e-sequence at positions si < S2 < ■ • • < Sk, the 
determinant is of a {k + 1) x (k + l) matrix with {i,j) entry {'^^~^^~^^^^^^) for < i,j < k 
with So = 0, Sk+i = n. 
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Example 1. If n = 8, ei = 1, £2 = £3 = €4 = 0, €5 = 1, eg = = 0, the matrix is 

/ 2 6 9 \ 

When 6=2, this is I 1 5 8 ) with determinant 9. Thus the chance of a carry at exactly 



.0 14, 

positions 1 and 5 when adding eight binary numbers is 9/2® = .03516. When b = 10, 
this chance becomes 1,042,470/10'^ = .0104. This is not so much smaller than in the 
binary case. 

Proof. (Of fact 5) We follow Stanley [67|, p. 69]. Let an{S) be the number of sequences 
with descent set contained in S and Pn{S) the number of sequences with descent set 
equal to S. Then 

aniS) = PniT), f3n{S) = ^ (- 1) l^"^' (T) . 

TCS TCS 

There is a simple formula for a„(S'): the number of weakly increasing sequences with 
entries in {0, 1, . . . , 6 — 1} of length / is i^^Si) using the usual stars and bars argument 
(place 6—1 bars into / + 6 — 1 places, put symbol to the left of the first bar, symbol 1 
between the first and second bar, and so on, with symbol 6 — 1 to the right of the last 
bar). Then 



an{S) 



si + 6-l\/s2-si + 6-l\ fn- Sk + h-l 
b-l )\ 6-1 j I 6-1 



because any such sequence of length n is constructed by concatenating nondecreasing 
sequences of length si, S2 — si, . . . , n — s^. Therefore 

/?n(S) = (-l)'~'/(0, h)f{h.i2) ■ . . f{i,, k + 1) 

l<il<i2<---<ij<k 

With f{i,j) = ('^ll+'"') (with SO = 0, Sk+1 = n). From Stanley's discussion, l3n{S) is 
the determinant of the (k+l) x (k+l) matrix with entry f{i,j+l), < i, j < k. □ 

Example 2. From the above development, the chance that the sum of n base 6 digits 
has no carries is ("|^7^)/6". When n = 8, 6 = 2, this is 9/2® = .03516. When 
n = 8, 6 = 10, this is ( g^) /lO® = .00024. We feel lucky when adding eight numbers with 
no carries. These calculations show that we are lucky indeed. 



3 Determinantal point processes 

Let A' be a finite set. A point process on A" is a probability measure P on the 21'^ 
subsets of X. For example, if A" = {1, 2, . . . , n — 1}, a point process is specified by 
recording where the carries occur in adding n base 6 numbers as in [Section "21 One 
simple way to specify P is via its so-called correlation functions. For A <^ JY, let 

p{A) = P{S : S^A}. 



4 



The collection of numbers {p{A)} uniquely determine P using inclusion/exclusion. A 
point process is determinantal with kernel K(x,y) if 



p(^) = det {K{x,y))^ y^^. 

On the right is the determinant of the \A\ x |^| matrix with {x,y) entry K{x,y) for 
x,y e A. 

Determinantal point processes were introduced by Macchi [50[] to model the distribu- 
tion of fermions. Se e [itI] for a textbook account, [63] for a survey, [4l| | for probabilistic 
developments and [4^] for many combinatorial examples. We warn the reader that 
these last references essentially deal with symmetric kernels K{x,y) whereas all of our 
examples involve non-symmetric kernels. 

There has been an explosive development of determinantal point processes over the 
past few years. Some reasons for this include: 

(a) A raft of natural processes (including all those in the present paper) turn out to 
be determinantal. These include: the eigenvalues of various ensembles of random 



matrices (Dyson [25|, |26|, |27[ , Mehta [55[ ) ; the presence or absence of edges if a 
random spanning tree is chosen in a iTaph (Burton-Pemantle 0, Lyons 0); 
random tilings and growth models (Johansson [3]); the structure of random par- 
titions chosen from a variety of measures (e.g., the Plancherel measure on the 
symmetric group Sn) (Borodin-Okounkov-Olshanski [H]); random dimers on bi- 
partite planar graphs (Kenyon [3]); and the zeros of a random analytic function 
E^=o«n^" With i.i.d. complex Gaussian coefficients a„ The Hst goes on 

extensively. 

(b) Specifying a kernel may be much easier than specifying a measure on all 2''^ 
subsets. Further, the correlation functions allow easy computation of quantities 
of interest. For example, if is 1 or as the random set includes x or not: 

£(XJ = A-(l,:lO, 

Co,(x,.x,) = de. (^^^ f^- f^ - K(...my,y) 

= -K{x,y)K{y,x). 

If K{x,y) = K{y,x), then the correlation is < 0, a distinctive feature of determi- 
nantal point processes. 

(c) If the matrix {K{x,y))x,y^x has all real eigenvalues {Xx}x&x-, many theorems 
become available for A^, the total number of points in a realization of the process. 
These include: 

• A is distributed as a sum Yl,x&x with {Y^} independent 0/1 random vari- 
ables having PiY^ = 1) = A^. 

• Let /i and a denote the mean and variance of A. These are available in terms 
of K from (b) above. Then the following refined form of the central limit 
theorem holds: 



27r 



oo 



.80 

< — 
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• The probabiUty density P{n) = P{N = n) is a Polya-frequency function. In 
particular the density is unimodal. Further, |mode — mean| < 1. 

• There is a useful algorithm for simulating from K (here assuming K{x, y) = 
K{y,x)). 

Most of these properties follow from the fact that the generating function of A^, 

E{x^) = +{x-l)K) 



has all real zeros. See Pitman 57[ or Hough, et al. [41] for details. 

(d) If Kn is a sequence of such kernels converging to a kernel then, under mild 
restrictions, the associated point processes converge. See Soshnikov [65| or Jo- 
hansson and Nordenstam [45[ for precise statements and remarkable examples. A 
simple example is given in [Section "51 (Example [TT]l . 

A main result of our paper is that the carries and descent processes, and indeed any 
one-dependent point processes on the integers, is determinantal. The special case of 
carries is stated here; a proof of the general case is in [Section 7l 

Theorem 3.1. The point process of carries and descents Pn given in \Section S\ is de- 
terminantal with correlation kernel K{i,j) = k{j — i) where 

y k(m)t"' = -V. 

^ ' i-(i-t)b 



Example 3. Take h - 
to k{rri)c!^ and K{i,j) 



2. Then 



Replacing t by ct changes k{m) 
^K{i,j)c^. Conjugating K hy a 



— 2t(l-t/2) ■ 

k{j — i) to c^~^k(j — i) = c~ 
diagonal matrix does not change determinants or the correlation functions. Thus setting 
4 = 2;, the generating function becomes 



\[- + l + z + z^ + 
4 V 2; 



When e.g., n — 1 = 5, the kernel is 



11111 
11111 
1111 
111 
1 1 



(3.1) 



From this, p{A) 



This is a manifestation of the 



HA has two consecutive entries 

(i)'^' otherwise 

one-dependence together with the fact that (for 6 = 2), two consecutive carries are 
impossible. It follows that the possible configurations are all binary strings of length 
n — 1 with no two consecutive I's. This is a standard coding of the Fibonacci number 
Fn+i (if 1, 1, 2, 3, 5 are the first five Fibonacci numbers). For example, when n = 4, the 
possible configurations are 000, 001, 010, 100, 101. 
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The eigenvalues of any finite matrix of the form (j3.ip are real. This follows from the 
fact that these matrices are totally positive, see {28], and the well known fact that the 
eigenvalues of a totally positive matrix are real (and nonnegative) , see e.g. Corollary 
6.6 of [1]. 



The classification of 28| also implies that the correlation kernels for 6 > 3 are not 
totally positive (for large enough matrix size), and we do not know if their eigenvalues 
are real. 

Example 4. Take 6 = 3. A straightforward expansion leads to 
^ 1 1 



k{m)f 



m=—oo 



l + t + ^ + ^- + ^- 



3M 1 - 1 + V 



3t 1 + 



27 

Thus 

• k{n) = for n < —1, 

• M-i) = i 

• k{n) = for n = 4 + 6j, < j < 00. For other n > 2, 

71+1 I "+3 I 

• k{n) = (-l)L-6-J (1) L 2 J 2'^("), where 6{n) = 1 if n = 1 mod 6 and else. 
For example, the first few values are: 
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It is amusing to see that the two-variable kernel K(x, y) = k{y — x) reproduces the 
correlation functions found for this problem in [Section "21 (Fact 4) . Thus for n = 7, 
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4 One-dependent processes 

Let {Xj},^Q be a binary stochastic process. Here < < 00. The process is one- 
dependent if {XjYjZ^Q, {Xj}j^^^_^i are independent for all i. A natural example: let 

{Yi}f^^ be independent uniform random variables on [0, 1]. Let h : [0, 1]^ {0, 1} be 
a measurable function. Then Xi = h(Yi,Yi^i), < i < A, is a one-dependent process 
called a two-block factor. Much of the literature on one-dependence is for stationary 
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processes. However, if Yi are independent but not identically distributed or h = hi 
depends on i, the associated Xi is still one-dependent. Our main results in [Section Tl 
are for general one-dependent processes. We first describe some results for the stationary 
case. 



Stationary one-dependent processes 

A book-length review of this subject by de Valk [t^ collects together many of the results, 
so we will be brief, focusing on later research. One of the first problems studied is the 
question of whether all one-dependent processes are two-block factors. Counterexamples 
were found [t!] though it appears that most "natural" one-dependent processes are 
two-block factors. For example, let / G ^^(tt) be a function on the unit circle. Let 
K{i,j) = f{i — j), i,j G Z, / the Fourier transform. Macchi [s^] and Soshnikov 
show that K{i,j) is the kernel for a determinantal point process on Z. Lyons and Steif 



[49| and Shirai and Takahashi [62, l63|] develop remarkable properties of the associated 
point process. If f{9) = ae~*^ + h + ce'^, the point process is one-dependent and 
shows it is a two-block factor. See also [13l |. 

A stationary one-dependent process is determined by the numbers ai = P{Xi = 
X2 = ■■■ = Xi_x = 1), ai = 1. Indeed, if / = /i, . . . , e = e^+i, . . . , e^-i are binary 
strings, P(/, 0, e) + 1, e) = P{f)P{e). From this P(/, 0, e) = P{f)P{e) - P{f, 1, e). 
By similar reductions, the probability of any pattern of occurrences can be reduced 
to a polynomial in the Oj. For example, P(0,0, 0) = 1 — 802 + a| + 203 — 04. There 
is a remarkable, simple formula for this polynomial as a minor of a Toeplitz matrix. 
We learned this from [58, Chap. 7]. A generalization to the non-stationary case is in 
[Section 7l 

Theorem 4.1. For a stationary one- dependent process, let Oj = P{Xi = X2 = ■ ■ ■ = 
= 1), ai = 1. Let tl, . . . ,tn-i be a binary string with j zeros at positions S = 
{si < • • • < Sk} ^ [n — 1]. Then 



P 



(tl, • • • , tn-i) = det (a^^.+i_s,)^=o ■ 



Here the determinant is of ak + 1 by k + 1 matrix, and one sets sq = 0, s^+i = n, oq = 1, 
and Oi = for i < 0. 

Proof. Theorem 14.11 follows from Proposition 2.2 and Theorem 4.2 in [s^ . Chap. 7], after 
elementary manipulations. □ 

Remark: It follows from the discussion on fssi, p. 140] that if all of the determinants 
(for all n and subsets S) in Theorem 14.11 are non-negative, then a stationary, one- 
dependent process with these Oj values exists. 

Example 5. For any n, the sequence with n — 1 ones has k = zeros. The relevant 
matrix is then 1 by 1 with entry a„,, giving this probability. 

Example 6. To compute P{0, 0, 0) one applies the theorem with n = 4 and S = 
{1,2,3}. It follows that 



P(0,0,0) = det 



1 02 as 04 

1 1 02 03 

1 1 02 

11 



9 

3a2 + 02 + 203 — 04, 
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as above. 



Our next result expresses the determinant in Theorem 14.11 as a skew-Schur function 
of ribbon (also called rim-hook) type. Background on skew-Schur functions can be found 
in 51, Sec. 1.5]. For further references and recent work, one can consult or 47]. We 
will use the elementary symmetric functions in countably many variables defined as 
eo = 1 and 

\ ^ 

l<ii<i2<-<ir 



From [5l|, p. 20] these are algebraically independent, so there is a homomorphism of the 
ring of symmetric functions to M which sends each to Oj from Theorem 14.11 



Theorem 4.2. With notation as in Theorem \4.1\ let A and fi be the partitions defined 
by 

Xi = n — Si-i — k + i — 1 , = n — Si — k + i — 1, 1 < i < k + 1. 

Let X' , fj,' denote the transpose partitions of X and fi and let Sy/f^' denote the correspond- 
ing skew-Schur function, obtained by specializing the elementary symmetric functions Ci 
to equal Oj. Then 

P{tl, ■ ■ ■ ,tn-l) = Sy/fj,'- 

Proof. First note that the Aj,;Uj defined above are all non-negative. From Macdonald 



sy/f,' = det(eA,-^^.-,+j)^+ii, (4-1) 



51|, p. 71], 

and Cr is the rth elementary symmetric function. When each is specialized to equal 
Or, the quantity sy/^i becomes 



det{ax,-f.^_i+j)^+l^ = det(a,^._,,_J^+i„ 
as desired. □ 

One reason why Theorem 14.21 is interesting is that there are non-obvious equalities 
between ribbon skew-Schur functions. The paper characterizes when two ribbon 
skew-Schur functions are equal; analogous results for more general skew-Schur functions 
are in [gQi]. In particular, combining the results of [3j] with Theorem 14.21 one immediately 
obtains the fact that a stationary one-dependent process is invariant under time reversal, 
i.e. 

P{ti, ■ ■ ■ ,tn) = P{tn, • • • ,ti) 



for all ti. For another proof of invariance under time reversal, see 58|, p. 139]. 

The determinantal formulae of Theorems 14.11 and 14.21 suggest that a determinantal 
point process is lurking nearby. This is indeed the case. In [Section 71 (see Corollary (73]) 
we prove the following. 



Corollary 4.3. A stationary one- dependent process as in Theorem \4-l\ is determinantal 
with kernel K{x, y) = k{y — x) with 

oo 

Y,Kn)z^ = -l/Y.<',zK 
nez i=i 
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A remarkable development, connecting stationary one-dependent processes to alge- 
bra appears in [ss!, Chap. 7]. They consider a graded algebra © © ^2 + • • • with 
Aq = k, a ground field, and AiAj C Ai^j. For example, the space k[x] of polynomials 
in one variable has Aq = k, and Ai spanned by x*. They assume that each Ai is a finite 
dimensional vector space of dimension dim(Aj). The algebra is quadratic if = k, 
the algebra is generated by elements of Ai, and the relations defining the algebra are 
in ^2- For example, the commutative polynomial ring k[xi, . . . ,Xn] is generated by 
xi, . . . , Xfi and the quadratic relations = 0. Note that algebras need not be 

commutative. 

A technical growth condition on the dim(Aj) which we will not explain here yields 
the Koszul algebras. These include many natural algebras occurring in mathematics 



(see 3j] for a survey). As a simple example, consider the commutative polynomial ring 
generated by xi,X2, ■ ■ ■ ,Xn with the additional relations XiXj = for £ E, with 
E the edge set of an undirected graph on {1,2,..., n| (loops allowed). This is Koszul 



|58| . Chap. 2, Cor. 4.3]. Pohshchuk and Positselski [58|, Chap. 7, Cor. 4.3] prove the 
following remarkable result. 

Theorem 4.4. To every Koszul algebra A one can assign a stationary, one-dependent 
process via 



Part of the reason that Theorem 14.41 is substantive is that the quantities P{Xi = 
X2 = ■ ■ ■ = Xi-i = 1) can not assume arbitrary values; the determinants in Theorem 
14.11 must be non-negative. 

Example 7. Consider the commutative polynomial algebra generated by xi, X2, ■ ■ ■ Xb 
and the additional relations = 0, 1 < i < b. The degree i part Ai is spanned by 
square free monomials and has dim(ylj) = (^), < i < 00. From Fact 2 of [Section 2^ 
we see that the associated one-dependent process is precisely the carries process of mod 
b addition from our introduction. Note that since only the dim(^j) matter, the algebra 
in this example can be replaced by the exterior algebra generated by xi,X2, ■ ■ ■ ,Xb and 
the relations XiXj = —XjXi for all i,j. 

Example 8. Consider the commutative polynomial algebra generated by xi,X2, ■ ■ ■ , x^. 
The degree i part Ai is spanned by monomials of degree i and has dim(Aj) = (^^^~^)- 
From Example [2] of [Section "2| the associated one-dependent process is precisely the 
complement of the carries process of mod b addition. More generally, in [Section "Tl 
we show that the particle-hole involution of a (perhaps non-stationary) one-dependent 
process with respect to any subset is one-dependent. 

Example 9. A PBW algebra is an algebra of the form 

k[xi,--- ,Xn]/{xiXj : {i,j) £ S C [l,nf), 
where the variables xi,--- ,Xn do not commute. Here it is useful to think of S as 



a directed graph with vertex set {h--- ,n}, loops allowed. Prom [58j, p. 84], PBW 



algebras are Koszul. As noted in [53, Chap. 7, Prop. 5.1], the one-dependent process 
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associated to a PBW algebra can be described as follows: pick Ui,U2,-'' i-i-d- in 
{!,••• ,n} and let Xi = h{Ui, f/j+i), where h{i,j) = if and only if G S. Indeed, 
the dimensions satisfy dim(74o) = l,dim(^i) = n, dim(Aj) =number of paths of length 
i — 1 in the complement of S, and so the chance of i — 1 consecutive I's in the point 
process is equal to dim(Aj)/ dim(^i)*. 

The above argument shows that PBW algebras give rise to two-block-factor pro- 



cesses. In fact they are dense in the set of all two-block factors 58|, Chap. 7, Prop. 5.2]. 
For a generalization of these PBW processes, as well as a simple argument that they 
are determinantal, see Remark 5 after Theorem 15.31 

ective space 



Example 10. Consider 2n points xi,X2, ■ ■ ■ , X2n in general position in pro 

mn 



The coordinate ring of this projective variety is known to be Koszul [58|, p. 42] with 
dim(^o) = l,dim(Ai) = n + l,dim(^i) = 2n for i > 2. Thus P{Xi = X2 = ■ ■ ■ = 
= 1) = 2n/(n + 1)* defines a one-dependent process. When n = 1, this is fair coin 
tossing. When n = 3, it is the carries process for multiplying quaternions introduced in 
[Section 6l 

For general n, these one-dependent processes can be described as follows (and illus- 
trates what one might consider to be a "probabilistic" description). Let {1,2, ••• ,n} 
be ordered cyclically and let * be another symbol. Choose Ui i.i.d. in {*, 1, • • • , n}. Let 
Xi = h(Ui,Ui-^^i) where h{*,x) = for all x, h{x,*) = 1 = h{i,i + 1) for x ^ *, and 
= otherwise. To see that this works, note that ^(-'^i = 1) = 2n/{n + 1)^, since 
one has to choose x different from * first and then * or x + 1 next. Similarly, for the 
chance of i — 1 I's in a row, the first choice can be anything but *, then the next i — 2 
choices are determined and the last can be one of two. 

We are certain that natural Koszul algebras will lead to natural point processes. We 
note further that [s^. Sec. 7.6] shows how natural operations on algebras preserve the 
Koszul property. These include the operations of union and complements that we work 
with in [Section 71 



5 Descents in permutations 



A permutation a G 5^, the symmetric group, has a descent at i if a{i) > a{i + 1). The 
set of such i forms the descent set D(a). If the base b in previous sections is large and 
n stays fixed, a string of n digits will have no repeated values and the descent theory of 
[Section "21 and [Section "31 above becomes descent theory for random permutations. This 
is a venerable subject. Stanle y [Gq . |68 ] reviews the basics. Some modern highlights 
are Solomon's descent algebra [6j], the connections with the free Lie algebra (30, '61|, 
Gessel's theory of enumerating permutations by descents and cycle structure [22, ,37i ], 
quasi-symmetric functions and the theory of riffle shuffling 35|, l70|] . Any Coxeter group 
has its own descent theory 6j|. There is some indication that arithmetic carries can be 
carried over [2o| . 

This section introduces three examples where the descent set can be shown to be 
a one-dependent determinantal point process: uniform choice of permutations, non- 
uniform choice from the Mallows model, and independent trials. Then it shows that the 
closely related notion of the connectivity set of a permutation also yields a determinantal 
process. Finally, connections with binomial posets are mentioned. 
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Example 11. The descent set of a uniformly chosen permutation. Consider 
the formula of ISection "2] Fact 5 for the chance that a random 6-ary string of length n 
has descents exactly at 5 C [n — 1]. Passing to the limit, as 6 ^ oo using (''^^Y^) ~ ^ 
gives a classical formula for the chance that a uniformly chosen a € Sn has descents 
exactly at S: 

Pn {D{a) = S) = det [1/ - s^)!] . (5.1) 

For S = 1 < si < S2 < ■ ■ ■ < Sis < n — I, the determinant is of a (/c + 1) x (/c + 1) matrix 
with entry l/(sj+i — Sj)!, (z,.? ) G [0, fc] x [0, fc] and sq = 0, s^+i = n. This formula 



is originally due to MacMahon [52] • See Stanley [67|, p. 69] or Gessel and Viennot [3. 
for modern proofs. 

Equation (|5.ip allows us to see that the descent set of a uniformly random permu- 
tation, treated as a point process, is determinantal. The proof follows from (|5.ip and 
Corollary 17.61 in ISection "7l 

Theorem 5.1. There exists a stationary point process on Z, call it P, such that its 
restriction to any interval of length n — 1 coincides with Pn of (|5.ip . The process P is 
one- dependent and determinantal. Its correlation kernel K{x,y) = k{y — x) with 



Remarks. 



1. Theorem 15.11 can also be seen by passing to the limit in Theorem 13. 11 Replacing t 
in the generating function 1/(1 — (1 — t)^) by —z/b and letting 6 oo gives (15. 2p . 



2. The Bernoulli numbers are defined by z/{e^ — 1) = Bnz"" /nl. These are 
very well-studied. It is known that i?2i+i = for i > 1, Bq = 1, Bi = —^,B2 = 
i, B^ = = ^) We see that k{m) = -Bm+i/{m + 1)!. 

3. From these calculations, the kernel K is 



K 



12 

1 

2 



J_ 

12 



— — 

n U 720 





The correlation functions of P (and P„) can be computed explicitly by passage to 
the limit from Fact 5 of lSection 21 Let A = U ^2 U • • • U be a decomposition 
of the finite set A C Z into disjoint blocks of adjacent integers. If \Ai\ = Oj, 



p{A) = n 



1 



i=i ^ * 



1)!' 



12 



5. For any n = 1,2, . . . , let Kn{i,j) = k{j — i) be the n x n top left corner block of 
K. Let d{TT) denote the number of descents of vr and let 



A 



n+l 



be the (n+l)st Eulerian polynomial (see e.g., [67|, p. 22]). Then item (c) preceding 
Theorem 13.11 yields that 

An+i{x) = det (/ + (x - l)Kn) (n + l)! 

It is known [3, p. 292] that An+i{x) has all real zeros a\,a2, ■ ■ ■ ,an- From the 
development in [Section 31 A^^^+i, the number of descents in a random permutation 
from Sn+i, is the sum of n independent Bernoulli random variables with success 
rates (1 — o-j)"^, 1 < i < n. One easily computes 



n+lj 



ti{K„ 



n 
2' 



Var(A^, 



n+l) 



tr (K, 

n 
2 



5n 1 
12 ~ 6 



n + 2 
12 



It follows by Harper's method ( [391] . [571] ) that one has a central limit theorem: 

Nn+l — f 



n+2 
12 



6. Chebikin [11[ studied the "alternating descent set" A{a) of a random permutation, 
defined as the set of positions at which the permutation has an alternating descent, 
that is an ordinary descent if i is odd, or an ascent if i is even. Combining Lemma 
2.3.1 of [lH with the argument used to prove (jS.ip . one obtains that for \S\ = k, 

Pn{A{a) =S) = det[Es^^,.sJis,+i - s.)!]u=o> 

where En is the nth Euler number defined by ^n>Q-^nZ^ /n\ = tan(z) + sec(2:), 
and So = 0, Sj^+i = n. Applying our Theorem 17.51 in [Section 7l it follows that 
that Pn is obtained by restricting to any interval of length n — 1 the stationary, 
one-dependent, determinantal process with correlation kernel K{x, y) = k{y — x) 
with 

^ k{m)z"' ^ 



1 — (tan(2;) + sec(z)) ^ 



7. Modern combinatorics suggests a host of potential generalizations of Theorem l5.ll 
Let P be a partial order on [n]. A linear extension of P is a permutation a £ Sn 
such that if i is less than j according to P, then a{i) < Let C{P) denote the 

set of linear extensions of P. For each a £ ^(P), let 



D{a) = {i : aii) > (T(i + 1)}, 

A(cr) = {i : a{i - 1) < cr(i) > a{i + 1)}, 



the descent set of a, 
the peak set of n. 



Choosing a uniformly in C{P) gives two point processes. This section has focused 
on ordinary descents, that is, P is the trivial poset with no restrictions. There 
are many indications that "natural" posets will give rise to determinant al p oint 
processes. For background and first results, see Brenti [3] or Stembridge 72]. 
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8. Finally, we mention that it would be interesting to find analogs of Theorem 15. II for 
other Coxeter groups. For the hyperoctahedral group Bn consisting of the 2"'n! 
signed permutations, descents can be defined using the linear ordering 

1 < 2 < • • • < n < -n < • • • < -2 < -1. 

Say that 

(a) a has a descent at position i (1 < i < n — 1) if a{i) > a{i + 1). 

(b) a has a descent at position n if <j{n) < 0. 



For example, —4 — 1 3 2 5 G -B3 has descent set {2, 3, 5}. Reiner [59[ shows that the 
chance that a random element of Bn has descent set {si, • • • , Sr} C {1, 2, • • • , n} 
is det(a(z, j))[j^o where 



Here sq = 0. Theorem 17.51 in [Section "71 shows that the resulting processes is 
determinantal (but not stationary). 

Example 12. Descents from non-uniform distributions on permutations. In 

a variety of applications, non-uniform distributions are used on permutations. A widely 
used model is 

Pe{a) = a G (5.3) 

Here, < ^ < 1 is a fixed concentration parameter, ctq G is a fixed location parameter, 
(i((T, (To) is a metric on Sn and Z = ^'^{'^''^o) is the normalizing constant. These 
are called "Mallows models through the metric d." When 9 = 1, Pi is the uniform 
distribution. For < ^ < 1, P51 is largest at ctq and falls off as a moves away from (Tq. 



See [15|, 23| or [5j] for background and references. 

Perhaps the most widely used metric is 

d{a, (To) = minimum number of pairwise adjacent transpositions 

to bring (T to iTq, (5.4) 
= / ((T(Tq"^) , the number of inversions in cra^^ . 

This is called "Kendall's tau" in the statistics literature. References, extensions and 



properties are in [19j, Sect. 4], [3l|], [53(]. We mention that the normalizing constant of 
the Mallows model through Kendall's tau is given by Z = YYi=i 0' P- 

With 6, (Tq and d fixed, one may ask any of the usual questions of applied probability 
and enumerative combinatorics: "Picking a permutation randomly from Pe{-), what is 
the distribution of the cycle structure, longest increasing subsequence, . . . ?" For general 
metrics, these are open research problems. Even the algorithmic task of sampling from 



Pe leads to difficult problems. See [2l|] or [It 

For Mallows model Pg through Kendall's tau (j5.4p and (To = id, we show that the 
descent set of it, D{a) = {i : a{i + 1) < cr{i)} forms a determinantal point process with 



simple properties. In what follows we use the g-analog notation that iq = ^^j- and 
n\q = nr=i ^9' t^^^ 6 = q as is conventional in combinatorial work in the subject. 
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Proposition 5.2. Let Pq (0 < q < I) be Mallows model through Kendall's tau ()5.3p . 
dED on Sn- 



(a) The chance that a random permutation chosen from Pq has descent set si < S2 < 
■ ■ ■ < Sk is 

1 



k 

det ' 



{Sj^ 



i,j=0 



with sq = and s^+i = n. 



(h) The point process associated to D{a) is stationary, one- dependent, and determi- 
nantal with kernel K(x,y) = k{y — x), where 



i-(e:^=o-'"M)" 



(c) The chance of finding k descents in a row is 2 ) j {k -\- Jn particular, the 
chance of a descent at any given location is q/{q + 1). The number of descents 
has mean ^{n, q) = ^^^^^^ and variance 



— — 9+^)" q +3q 1 ^ Normalized by 



its mean and variance, the number of descents has a limiting standard normal 
distribution. 



Proof. Part (a) of Proposition [5^2] follows from Stanley [671 , Ex. 2.2.5]. Part (b) follows 
from Corollary 17.61 in ISection 7l and elementary calculations. To compute the chance of 
k consecutive descents, note by stationarity that this is the chance of the permutation 

o-(i) = A;+2—z in S'fc+i under the measure on This probability is g( 2 )/(fc+l)!q, 

as needed. One calculates from part (b) that 

k{-l) = 1, k{0) = k{l) 



+ 1' ' ' (g + l)2(g2 + g + i)- 
Recall from the remarks preceding Theorem 13.11 that 

E{x^) = det (/ + (x - l)i^) , 
where N is the total number of particles. Thus 

E{\D{a)\) = tr{Kn-i) = (n - 1)A:(0), 

Var{\D{a)\) = tr{Kn-i - K^^i) = (n - 1)A;(0) - (n - l)k{Of - 2{n - 2)k{l)k{-l), 
and the proof of part (c) is complete. □ 
Remarks. 

1. From part (c), the correlation functions are given explicitly, just as in remark 4 
following Theorem 15.11 

2. It would be interesting to have a direct combinatorial proof of the one-dependence 
in Proposition! 
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Example 13. Descents for independent trials The carries in tiie carries process 
have an equivalent description in terms of descents in an independent and uniform 
sequence. In this example we generalize to descents in a sequence Yi,Y2, . . . ,Yn with 
Yi independent and identically distributed with P{Yi = j) = pj. Here < j < 6 — 1 
to match previous sections (and of course < Pj < IjSPj = l)i but in fact b = oo 
can be allowed. The descents are related to the order statistics for discrete random 
variables [3, [i^ but we have not seen this problem previously treated in the probability 
or statistics literature. Of course, if Yi are independent and identically distributed from 
a continuous distribution, the descent theory is the same as the descent theory for a 
random permutation. 

It is useful to have the complete homogeneous symmetric functions 

(-^l ) X2 , • • • , S^m) — ^ ^ Xi-^Xi2 . . . Xjj. . 

l<ii<---<ij.<m 

Theorem 5.3. IfYi,Y2, . . . ,Yn are independent with common probability P{Yi = j) = 
Pj, < j < 6 — 1, the descents in Yi, Y2, . . . ,Yn form a stationary one-dependent point 
process on [n — 1] with 

Pn{S) = det (^+i-s,)o<i.i<fe ^^^^ '^0 = ^' '^'=+1 = 
for S = {si < S2 < ■ ■ ■ < Sk}. The hj are evaluated at pQ, . . . ,pb^i. Moreover, this 



process is determinantal with K{i,j) = k{j — i) and Xli^-oo ^(O-^' 



1 



i-nto(i-p^^)' 



Proof. Let a„(S') = P„(des(yi, I2, • • • , ^n) ^ S) = hs^ hg^-s^ ■ ■ ■ Ksk- The determi- 
nant formula follows from the argument for Fact 5 in [Section "21 Now, Corollary 17.61 
implies that k{z) = , with e{z) = Ya^q^iz^ = IltdCl - Piz)~^ from e.g., 

Macdonald [ill, p. 21]. □ 



Remarks. 

1. Setting Pi = 1/6, k{z) = ■ This is equivalent to our result in Theorem 
13. H by changing variable z to z/b — > t, which doesn't change correlations. 

2. Prom the theorem, the associated point process is one-dependent. If Xj = 1 or 
as there is a descent at i or not, then clearly 

P{Xi = X2 = ■ ■ ■ = Xi = 1) = ei+i(po, • • • ,Pb-i), 

where er{xi,--- ,Xm) = Z]i<n<-..<v<m ^n^j2 ' ' ' Xi^ is the rth elementary sym- 
metric function. For example, 



i 

p{x. = x,^, = i) = 1-\J2p' + IEpI 

i i 

This allows a straightforward formulation (and proof) of the central limit theorem 
for the number of descents. As in [Section 21 simple expressions for the correlation 
functions are available in terms of P{Xi = X2 = ■ ■ ■ = Xi = 1). 
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3. As we show in Corollary 17.71 one can express PniS) as a skew-Schur function of 
ribbon type. For special choices of the pj's, tjiis symmetric function interpretation 
may facilitate computation. For example, |51|, Sect. 1.5] gives Sx/f^ = J^u^fiu^i^ 
with c^^ the Littlewood-Richardson coefficients and the Schur function. If 
Pi = zq^ (i = 0, 1, • • • ,b — 1, z a normalizing constant), then Sy is non-0 only if v 
has length < b, in which case 



^b+c(x) 

^ 11 1 _ fjh(x) 



see Macdonald [sT, p. 44] for an explanation of notation. 

4. Coming back to carries, if a column of numbers is chosen independently on 
{0, 1, . . . , 6 — 1} from common distribution {PjI^Zq, the remainder column after 
addition evolves as a Markov chain with starting distribution {pi} and transition 
matrix P(i,j) = pj-i (indices mod b). This is a circulant with known eigenvalues 
and eigenvectors. We do not know the distribution of carries. 

5. There is a sweeping generalization of this example which leads to a large col- 
lection of determinantal point processes. Let R be an arbitrary subset of [N] x 
[N]. Fix a probability distribution 9i,92, ■ ■ ■ ,0j\[ on [N]. A point process P„ 
on [n — 1] arises as follows: pick Yi,Y2, . . . ,Yn independently from {Oi}. Let 

Xi = ^ ^, l<i <n-l. Let S = {i : Xi = 1}. Using essentially 

I otherwise 

the arguments above, it is proved in Brenti Th. 7.32] that 

r R i'' 

Pn{S) = det h for 5" = {si < • • • < Sk}, sq = 0, Sk+i = n, 

L J i,j=0 

where hf = hf{9i,. ..,0n) ■= P0)- Here as usual, h(} = l,hf = for j < 0. 
This falls into the domain of Theorem 17.51 and its corollary. It follows that P„ is 
determinantal with kernel K{i,j) = k{j - i), ^(0^' = \«zj)~') ' 

f 1 if i < j 

Our examples are the special case R(i, j) = { ~ ■ We are certain that 

^ ^ V )j; 1 Q otherwise 



there are many other interesting cases. 



Example 14. Connectivity set of a permutation. 



Stanley 7l| studied a "dual" notion to the descent set of a permutation, which he 
called the connectivity set. For o" G S„, the connectivity set C{a) is defined as the set 
oi i, 1 < i < n — 1, such that a S x S'i+i^...^„. For example, the permutation 

(T = 312546 satisfies C{a) = {3, 5}. The connectivity set also arises in the analysis 



of quicksort, where it is called the set of splitters of a permutation 7^, Sec. 2.2]. 

We will prove that the connectivity set of a random element of Sn is a determinantal 
point process and determine its correlation kernel. For this the following two facts are 
helpful: 

1. A permutation a G Sn is called indecomposable or connected if C{a) = 0. Comtet 
[ij, Exer. VII. 16] shows that the number /(n) of connected permutations in Sn 
satisfies 

n>l 2^n>0^-^ 
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(An asymptotic expansion of f{n) is given in Exer. VII. 17]). 
2. [7l|, Prop. 1.1] Letting S = {si < S2 < ■ ■ ■ < Sfc} be a subset of {1, 2, • • • ,n — 1}, 
\{aeSn:SC Cia)}\ = si!(s2 - si)! • • • (sfc - Sfc_i)!(n - s^)!. (5.6) 

The following is our main result. 

Theorem 5.4. Let C{cr) denote the connectivity set of a permutation a chosen uni- 
formly at random from Sn- The point process corresponding to C{a) U {0, n} is deter- 
minantal, with state space {0, • • • , n} and correlation kernel K(x,y) satisfying 

KiO,y) = l ally 
K{n,y) = 5n,y ally 
K{x,y) = -^ 0,ny^x<y 

K{x, y) = jlr- -pi^ 0, n / X > y 

\x / \n-~x/ 

Note that in the statement of the theorem, 0, n are always points of the process. 

Proof. The first step is to observe that C{a) U {0, n} can be obtained as a trajectory of 
a certain Markov chain, started at 0, with transition probabilities 

[n — I)'. 

Here we take /(O) = 0, so P{i,i) = for all i. To prove (j5.7p . note that C{a) = 
{ii, • • • , ik} if and only if the following events Ei, - ■ ■ , E^+i occur: 

• Ei: {vr(l), • • • , 7r{ii)} = {1, • • • , ii} and tt restricted to {1, • • • , ii} is indecompos- 
able. 

• E2: {vr(ii + 1), • • • , vr(i2)} = {ii + 1, ■ ■ ■ , ^2} and vr restricted to {ii + 1, • • • , 12} is 
indecomposable. 



• Ek+i'. {vr(ifc + 1), • • • , vr(n)} = {ik + 1, • • • , n} and vr restricted to {ife + 1, • • • , n} 
is indecomposable. 

Letting f{n) denote the number of indecomposable permutations in Sn, the probability 
of El is clearly /(lll^pill. The probability of E2 given Ei is and the 

probability of £"3 given Ei,E2 is '^^^^^^i'^fj^^p^^' , etc., as claimed. 

Now 0, Thm. 1.1] implies that the point process on C{a) U {0,re} is determinantal 
with correlation kernel 

K{x,y) = 5o,, + Q{0,x)-Qiy,x) (5.8) 

where 

Q = P + p2 _^ p3 _^ . . . . 
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To compute Q, let denote the coefficient of in a series g{x), and note that 



I 

{n-l)lf{l-i) {n-j)\f{j-l) 
^ {n-iy. (n-/)! 

[n-iy. y l^n>on-^ J 

The last equality used (|5.5p . A similar computation shows that 



and thus 



{n-iy. \ E„>o^'^" 



(n-i)! 



>n>l 



Thus Q{i,j) = if ^ < J, and is otherwise. The theorem follows from this and 

dSS]). □ 

Remarks. 

1. The connectivity set C(cr) is a simple example of a determinantal process which 
is not one- dependent. Indeed, from P{1 G C{a)) = P{3 G C{a)) = 

and P(l,3 G = ^^^^ / P(l G C(a))P(3 G (7(a)). 

2. Unlike the other point processes considered in this paper, the expected number of 
points tends to as n — > cxd. Indeed, applying (|5.6p gives that 

n— 1 n—1 



E\C{ay = J2P{i€C{a)) = Y,-^^0 

i=l i=l \i) 



as n —>■ oo. 



Binomial posets 



In Doubilet-Rota-Stanley [2j|, binomial posets were introduced as a unifying mecha- 
nism to "explain" the many forms of generating functions that appear in enumerative 
combinatorics. Briefly, is a binomial poset if for every interval [x, y], the length of all 
maximal chains is the same (say n^y) and the number of maximal chains in an interval 
of length n does not depend on the particular interval. This number B{n) is called 
the factorial function. For the usual Boolean algebra of subsets, B{n) = n\. For the 
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subspaces of a vector space over a finite field, B{n) = n\q (g- factorial) . There are many 
further exarnples. 

Stanley [66j| has shown that many enumerative formulae generalize neatly to the 
setting of binomial posets. Some of these developments give new determinantal point 
processes. Here is an example. 

Example 15. (union of descent sets) 

Theorem 5.5. Pick cJi, C72, . . . , cr„ independently from Mallows model through Kendall's 
tau ()5.3p on the symmetric group Sn- Let S = {si < ■ ■ ■ < Sk} be the union of their 
descent sets. Then 

Pn{S) = det ^ ^ 



_{{Sj + l - SiY-qY 

where sq = 0, Sfc+i = n. The associated point process is stationary, one- dependent, and 
determinantal with K{x,y) = k{y — x) for k{z) = where e{z) = Yl^o^^ li^^-iY ■ 

Proof. The formula for Pn{S) is Corollary 3.2 of [66], so the theorem follows from one 
of our main results proved in [Section 71 (Corollarv l7.6p . □ 

When g = 1 and r = 2, e{z) = Iq{2^/z), the classical modified Bessel function. 



Feller [30|, §11.7] develops a host of connections with stochastic processes. 

Remark. The intersection of two independent point processes with correlation functions 
p^{A)^p'^{A) is a point process with correlation function p^{A)p'^{A). If both processes 
are determinantal, this is a product of determinants, but in general the intersection 
or union of determinantal processes is not determinantal. In Subsection 17. H we show 
that all one-dependent processes are determinantal. Since the intersection of descent 
sets of independent permutations is one-dependent, and taking complements preserves 
one-dependence, this gives another proof that the union process in Theorem 15.51 is 
determinantal. 

We have not pursued other examples but again believe there is much else to be 
discovered. 



6 More general carries 

Consider the quaternions Q% = {±1, iti, itj, itA;}. The center Z of Qs is {il}- Choose 
coset representatives X = A;} for Z in so any element can be uniquely repre- 

sented as g = zx. The coset representatives are multiplied by the familiar rule k r\ i r\ 
j r\ k so ij = k and kj = —i, etc. If we multiply gig2 ■ . . gk = {ziXi){z2X2) . . . (zkXk), 

the Zi £ {±1} can all be moved to the left and we must multiply xiX2 x^, keeping 

track of the "carries", here ±1. Evidently, if {gi} are chosen uniformly at random in 
Qs, both {zi} and {xj} are independent and uniform in Z and X. Thus we have the 
following problem: choose Xi,X2, . . . ,Xi^ uniformly in X and multiply as Xi, X1X2, 
{XiX2)X3, . . . , 

X2 
X3 

Xk 
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This gives a process of remainders and carries as in [Section Tl 



gives (— 1)^A; 



Example 16. 

k ■ k 

k 1 

i ■ i 

i 1 

k ■ k 

j ■ i 
k ■ j 

i k 

It is almost obvious that the carries form a stationary, one-dependent, two-block process 
with P{i — 1 ones in a row) = ji, 2 < i < oo. Further, the remainders in the second 
column are independent and uniform on {l,i,j,k} and there is a simple "descent" rule 
which determines the joint law of the dots (Example 1181 below). 

We also mention, by comparison with Example [10] in [Section "il that the carries 
process for the quaternions is the same point process that arises from the coordinate 
ring of 6 generic points in projective space P'^. 

One natural generalization where all goes through is to consider a finite group G 
and a normal subgroup contained in the center of the group. The factor group 
F = G/N has elements labeled 1,0", r, ... . We may choose coset representatives t{l) = 
1, t((T), t(r), . . . and any g = nt{a), uniquely. While sometimes t(o")t(r) = t{aT), in 
many cases this fails; but we may choose correction factors f{a,T) in N (often called 
a "factor set") so that t{a)t{T) = t{aT)f{a,T). Once t{a) are chosen, the f{o-,T) are 
forced. 

Example 17. If G = Cioo, N = Cio (thought of as a subgroup {0, 10, 20, . . . , 90}), 
the natural choice of coset representatives for G/N = {0, 1, . . . , 9} is t{i) = i ^ G. Of 
course, 



t{i)+t{j)=t{i+3) + f{i,j), with/(i,j) 



1 if i + i > 10 

if < « i < 10. 



It is natural to ask if the choice of cosets matters. To see that it can, consider Cio in 
Cioo and choose coset representatives as 0,11,22,43,44,45,46,47,48,49. The sum of 
11 and 22 requires a carry of 90. 

A lovely exposition of carries as co cycles is in ^4^. 

If 5 = nt{a) is chosen uniformly at random, then n and t{a) are independent and 
uniform. Multiplying a sequence of gi = riiti can be done by first multiplying the tj, 
keeping track of the carries, and then multiplying the rij and carries in any order. Of 
course, here the carries are in A^, not necessarily binary. 

Given ti,t2, ■ ■ ■ ,tk, we may form a two-column array with tj in the first column, 
the successive remainders ri, r2, ... in the second column and "carries" /i, /2, • • • fk-i 
(elements in A^) placed in between. 

Lemma 6.1. // coset representatives ti, 1 <i <k, are chosen uniformly, then 
1. the remainder process ri, 1 < i < k, is uniform and independent, and 
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2. the carries fi, 1 < i < k — 1, form a stationary, one- dependent, two-block factor. 
Proof. 

1. Since ri = ti, ri is uniform. Successive rj are formed by multiplying rj_i by 
ti. There is a unique choice of ti giving r^, so rj is uniform and independent of 
ri, . . . ,rj_i. 

2. Consider two successive remainders and the unique t giving rise to them: 

Ti-l 
it fi 

Since ti is uniquely determined, rj_itj = rif{ri^i,ti) is uniquely determined. It 
follows that the fi process is a two-block process: generate ri,r2... uniformly 
and independently, set fi = h{ri,ri^i), with h{r,r') = f{r,t) where r~^r' G 
tN determines t uniquely. Because two-block processes are one-dependent, this 
completes the proof. □ 

Corollary 6.2. With the notation of the lemma, define a binary process Bi, . . . , B]^_i 
as 

^ ^ fl ^//./id 
Then {Sj} is a stationary, one- dependent, two-block factor. 

Example 18. As in Example [M let G = Qs, N = {±1}, F = G/N ^ C2 x C2. 
The natural choice of coset representatives fc} gives rise to the following two- 

block representation: choose Ui, U2, . . . uniformly and independently in {1, i,j, k}. Let 
Bi = h{Ui,Ui+i) with h{l,x) = 1 (all x), h{x,l) = -l{x ^ 1), h{i,j) = h{j,k) = 
h{k,i) = —1, h{i,k) = h{k,j) = h{j,i) = 1, h{x,x) = 1, (x 7^ 1). 

Example 19. (dihedral group) Let Dq = {x,y : x'^ = y'^ = (xy)^ = !)• z = xy, this 
8-element group has center = = —1}. Choosing coset representatives l,x,y,z, 

the cosets multiply as 





1 


X 


y 


z 


1 


1 


X 


y 


z 


X 


X 


1 


z 


y 


y 


y 


—z 


1 


— X 


z 


z 


-y 


X 


-1 



From this, elementary manipulations show that P{i ones in a row) = 1/4*. Thus, the 
carries process is independent with P{Bi = 1) = 1/4, P{Bi = 0) = 3/4 for all i. 

We mention that Dg can also be represented as the extension of the normal subgroup 
C4 by the factor group C2. Here, coset representatives for C4 can be chosen so that 
there are no carries (i.e. the extension "splits"). 

Example 20. (extensions of C2 by Cm) A central extension of A^ = C2 by Cm is abelian. 
It follows that when m is odd, G = C2m is the only central extension of A^ = C2 by Cm, 
and when m is even, there are two central extensions C2 x Cm and C2m- The extension 
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C2 X Cm splits and choosing {(0, i), < i < m — 1} as the coset representatives, there are 
no carries. For C2m., with C2 = {0, m}, choose coset representatives 0, 1, 2, . . . , m — 1. 
Thus 

{1 \i i + j >m 
if < i + j < m. 



As for usual addition, there is a carry if and only if there is a descent in the remainder 
column. Thus 

n 



P{Bi = B2 = --- = 



Example 21. Let G = C2 x C2 x C2 and let N = {(0, 0, 0), (1, 1, 1)}. With coset 
representatives (0, 0, 0), (1, 0, 0), (0, 1,0), (1, 1, 0), there are never carries, but with coset 
representatives (0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), there are carries. This example shows 
that the one-dependent process i?i,i?2, • • • depends not only on G and N ^ but also on 
the choice of coset representatives. 

We have not embarked on a systematic study of carries for finite groups and believe 
that there is much more to do. We do note that by a result in the next section (Theorem 
I7.ip . the above processes, being one-dependent, are determinantal. 

The basic carries argument works for infinite groups as well. Let G be a locally 
compact group and H a closed normal subgroup. Suppose that H is in the center of 
G and that G/H is compact. Choose coset representatives t{(T) G G for o" G G/H. 
As in the finite case, these define factor sets f{(7,T) by t{a)t{T) = t{aT)f{a,T). Write 
g = nt. Since G/H is compact, it has an invariant probability measure. Choosing 
ti,t2, - ■ ■ independently from this measure and multiplying as above gives remainders 
and a carries process. Just as above, the remainders are independent and uniformly 
distributed in G/H and the carries process (with values in A^) is a one-dependent, 
two-block factor process. 

Example 22. A lovely instance of this set-up explains a classical identity. We begin 
with the motivation and then translate. Let cr be a permutation with number of descents 
d{a). It is known that, for a chosen from the uniform distribution on Sn, 

P {d{a) = j) = P (j < [/i + ... + [/„< i + 1) . (6.1) 

On the right, Ui,U2, - ■ ■ ,Un are independent uniforms on [O^JJ. The density and distri- 
bution function for [7i + - • • + [/„ was derived by Laplace; see [30] • Foata [sl] proved (16.1 ) 
by combining Laplace's calculation with an identity of Worpitzky. Richard Stanley [69[] 
gives a bijective proof involving an elegant dissection of the n-dimensional hypercube. 
Jim Pitman [s^] gives the following "proof from the book" which is a continuous version 
of our carries argument from [Section It form two columns 

Ui Vi 
U2 V2 

On the left are independent uniforms on [0, 1]. On the right are their remainders when 
added mod 1; so Vi = Ui, V2 = Ui + U2 (mod 1), • • • . The Vi are similarly independent 
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uniforms on [0, 1]. Place a dot at position i every time the partial sum Ui + ■ ■ ■ + f/j+i 
crosses an integer. Call these dots carries. As in the discrete case, there is a dot at 
position i if and only if there is a descent V^+i < Vi. The number of dots is the integer 
part of Ui + ■ ■ ■ + Un and also the number of descents, proving (j6.ip . Of course, the 
distribution of the descent process is the same as the carries process. 

In the language of group theory, let G = M, = Z, and G/N = Si, the circle group. 
Choose coset representatives as [0, 1) and factor sets in {0, 1}. 



7 Proofs and generalizations 



This section proves two of our main results for general one-dependent processes (we do 
not assume stationarity in this section). In Subsection 17.11 it is shown that all one- 
dependent point processes on Z are determinantal, a result which is new even in the 
stationary case. Subsection 17.21 proves that a point process P on a finite set X, with P 
given as a certain-shaped determinant, is one-dependent and determinantal. This covers 
quite a few examples from previous sections and is particularly useful in situations (such 
as Example [T2|) where the one-dependence is not apriori obvious. 



7.1 One dependent processes are determinantal 

For a random point process on a discrete set X, we define the correlation function p by 

p{A) = P{S : S^A}. 

Then one dependence on (a segment of) Z is equivalent to the condition that p{X[JY) = 
p{X)p{Y) whenever dist{X, Y) > 2. 

Theorem 7.1. Any one-dependent point process on (a segment ojTL) is determinantal. 
Its correlation kernel can be written in the form K{x,y) = 

ifx-y>2 
—1 if X — y = 1 

E'rT\-^Y-' E.=/o<h<...<^.=.+l P(i^O, h)) p{[hM))--- P{[lr-l,lr)) ^f X < y 

Here the notation [a, b) stands for {a, a + 1, • • • ,b — 1}. 

For example, K{x, x) = p{{x}), K{x, x + 1) = p{{x, x + 1}) — p{{x})p{{x + 1}), etc. 
Proof. By one-dependence, it is enough to verify that with K as above, 

det[K{x + i,x + j)\f]Zo = p{[x,y + l)) 

for any x < y. 

We use induction on y — x. For y = x the statement is trivial. Otherwise, one has 
that 



det[K{x + i,x + j)]l = det 



K{x,x) K{x,x + 1) ... K{x,y) 
-1 K{x + l,x + l) ... : 

-1 ■•. : 

-1 K{y,y)_ 
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When expanding this determinant, various numbers of — I's from the subdiagonal 
can be used. Observe that if we do not use the —1 in position — 1), then we 
can compute the corresponding contribution, because if we replace that —1 by 0, the 
determinant sphts into the product of two each of which is computable by the induction 
hypothesis. 

Similarly, if we insist on not using several —I's, then the contribution is the product 
of several determinants. Thus we obtain by inclusion-exclusion that 



det[K{x + i,x + j)]l-'' = K{x,y)+ ^ p {[x,l)) p i[l,y + 1)) 

x<l<y+l 

- Yl P{[x,li))p{[h,l2))p{[l2,y + i)) + 

x<li<l2<y+l 



(7.1) 



The K(x, y) term corresponds to using all —I's, the sum over I corresponds to not using 
a —1 at least at location (/,/ — 1), the sum over li,l2 corresponds to not using a — 1 
at least at locations — = 1,2, etc. Our definition of K(x,y) is such that the 

right hand side of (j7.ip is equal to p{[x,y + 1)), as desired. □ 

To state some corollaries of Theorem 17.11 we use the concept of particle-hole invo- 
lution. Essentially, given a subset M oi X, the involution maps a point configuration 
S C X to S A M (here A is the symbol for symmetric difference). This map leaves 
intact the particles of S outside of Af, and inside M it loses the particles and picks up 
the "holes" (points of J\f free of particles) . 

Corollary 7.2. The class of one-dependent processes is closed under the operations of 
particle-hole involution on any fixed subset, intersections of independent processes, and 
unions of independent processes. 

Proof. For intersections, the claim follows from the definitions. 

For particle-hole involutions, note that the property of being one-dependent follows 
from the kernel begin on the second subdiagonal and below. On the other hand, for 
determinantal point processes the "complementation principle " 0, Sec. A.3] says that 
the particle-hole involution can be implemented by the following change in the kernel: 



'A 


B 




'A B ' 


C 


D 




-C I-D 



where the block structure corresponds to the splitting into the noninverted and inverted 
parts. Clearly, this keeps the property of having O's on the second subdiagonal and below 
intact. 

Finally, unions can be reduced to intersections by the particle-hole involution. □ 

Remark. The class of determinantal point processes is not closed under intesec- 
tions/unions. 

Let us now look at the translation invariant case. Then p ([x, y)) = Py-x is a function 
oi y — X only, and pk is the chance of k consecutive ones. Set 

R{z) = l + z + Y,PkZ^^^. 

k>l 
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Corollary 7.3. In the translation invariant case, the kernel K{x, y) = k{y — x) is also 
translation invariant, and 



Proof. 



R{z) 



1 - R{z) 



Z 1 + PIZ + P2Z'^ + 

1 - z'™^" +\Y.p- 



1 1 

-- + - 

z z 



m>l 



E 

m>l 



.m>l 



E 

, m>l 



PrnZ 



+ 



which agrees with the formula of Theorem 17.11 



□ 



The generating function R{z) also behaves well with respect to the particle-hole 
involution on Z. 

Proposition 7.4. The particle-hole involution on Z with N = TL replaces R{z) by 
l/R{-z). 



Proof. We have 



1 



R{-z) 



1 



l-l/R{-z) R{-z)-l l-R{-z)' 
Hence by Corollary (731 changing R{z) to R{z) = l/R{—z) leads to the following change 



in the correlation kernel: 



k(n) = 5q., 



-lTk{n). 



This is equivalent to K{x,y) = 6x,y — \_]\y K{x, y). This is the same as K ^ I — K, 



which corresponds to the particle-hole involution [ai, Sec. A.3]. 



□ 



Example 23. The descent process on Z (Example II ip corresponds to pm = (rn+i)\ ' 
thus R{z) = e^. The particle- hole involution is given by R{z) = l/e~^ = e^, which is 
the same. 

Example 24. The intersection of r independent descent processes on Z corresponds to 
Pm = (^^j^rjjTF) hence 



m>0 



{m\y 



The correlation kernel is given by "^nez k{n)z'^ 



Example 25. The union of r independent descent processes is the particle-hole invo- 
lution of Example [Ml Thus Rjj\z) = l/Rj^\-z) and 

1 



Yk{n)z^ 



1 - 1/R, 



Replacing z by —z doesn't affect correlations, and we recover Example [TCI 
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7.2 A class of determinantal processes 

For any n = 2, 3, . . . , consider a probability measure Pn on all subsets S = {si < S2 < 
■ • • < Sfc} ^ [ra - 1] given by 

PniS) = h{n) det [e{si, J^.^^ , (7.2) 

for some /i : N ^ C and e : N x N ^ C with the notation sq = 0, Sk+i = n. We assume 
that e{i,j) = for i > j, and that e{i,i) = l,e{i,i + 1) > for all i. 

Theorem 7.5. If (|7.2p holds for some fixed n with e{i,j) = for i > j and e{i,i) = 
l,e{i,i + 1) > for all i, then Pn is a determinantal, one-dependent process with 
correlation functions 

p{A) = Pn{S : S ^ A} = det [K{ai, aj)]'^.^-^ for A = {ai, a2, . . . , Um}, 

with correlation kernel 

K{x,y) = 6^,y + {E~^)^,y+i, (7.3) 
where E is the upper triangular matrix E = [e{i — l,j)]",=i, 



E 



e(0, 1) 



e(0,2) 
e(l,2) 



e(0,n) 
e(l,n) 

e(n — 1, n) 



In addition, h{n) = (det E)-^ = (e(0, l)e(l, 2) . . . e(n - 1, n))-^ . 

Corollary 7.6. Assume further that e{i,j) = e{j — i). Then the point process is sta- 
tionary. If e{z) = X^/^o^(0-^') then K{x,y) = k{y — x) and 

k{z) := k{l)z^ = 



Proof of Theorem \ 7. 5| Set L = [e{i — 1, j) + -^j^. Thus L appears as 



L 



e(0,l) e(0,2) 
1 e(l,2) 
1 



e(0,n) 



1 e(n — 1, n) 
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For any function / : {1, . . . , n — 1} ^ C, by (|7.2p . 



E I n /(^« 



Pn ({Si < S2 < • • • < Sfc}) /(Si) . . . /(Sfc) 
+ 1, . . . ,Sfe + 1 



0<si<---<Sfe<n 
fc=0,...,n— 1 



det 



0<si<-<Sfe<n 
fc=0,...,n-l 



L 



si,...,Sfc,n 



/i(n) • det 



7(l)e(0,l) /(2)e(0,2) 
/(1)-1 /(2)e(l,2) 

/(2) - 1 



/(n- l)e(0,n - 1) e(0,n) 
/(n-l)e(l,n-l) e(l,n) 

/(n-l)e(2,n-l) e(2,n) 



/(n — 1) — 1 e(n — 1, n) 



/i(n) • det 



-1 



L 



7(1) 



fin -I) 



1 



This is the generating functional of P„. In the second equahty, the determinant is of 
the minor of L with rows 1, si + 1, • • • , Sfc + 1 and columns si, • • • , Sk,n. In the third 
inequality, the 2"~^ possible summands correspond to choosing which of the first n — 1 
matrix columns use the —1 coming from f{i) — 1 in the determinant expansion. 

If / = 1 + 5, the generating functional can be expressed in terms of the correlation 
functions: 

1+ X] Pm{si,...,Sm)gisi)...g{Sm) 

Sl<---<Sm 

m=l,2,... 



E 



SieS 

( 



h(n) det 



5(1) 



gin - 1) 



h{n) det(£^) • det 



I + E-'^L- 



\ 



5(1) 



g{n - 1) 



This holds for all g. First take g = to see h{n) det(ii^) = 1. Next note the expansion: 
if M is n X n, then det[/ + M] = X]5det(M(5)) with the sum over all 2" subsets of 
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[n], and M{S) the minor with rows and columns in S. This gives that Pm{si, ■ ■ ■ ,Sr, 
det[K{si,Sj)]'^j^-^, where K{x,y) = {E-^L)^y. Then 



E-^L = E 



-1 



V 



1 



+ E 



I + E 



-1 



1 



□ 



and the proof of (j7.3p is complete. 

Finally, we note that since K{x, y) vanishes below the first subdiagonal, P„ is one- 
dependent. 

Proof of Corollary. Here E is a Toeplitz matrix with symbol (e(z) — l)/z. It is triangu- 
lar, so E~'^ is a Toeplitz matrix with symbol z/{e{z) — 1). Thus K is a Toeplitz matrix 



with symbol k{z) 



1 + 



f(IY=T l-ieiz)) 



□ 



We also note that in the translation invariant case, there is an expression for Pn{S) 
in terms of skew Schur functions of ribbon type. 



Corollary 7.7. As in Corollary\7. 
fi be the partitions defined by 



assume further that e{i,j) = e{j — i). Let A and 



A,; 



n 



Si-l 



k + i - I , /ij 



n 



k + i -I, I <i <k + 1. 



Let X' , fj,' denote the transpose partitions of X and fi and let sxij^i denote the correspond- 
ing skew-Schur function, obtained by specializing the elementary symmetric functions ei 
to equal e{i). Then Pn{S) = ^^yj^sa'/t'- 

Proof. By assumption, formula (|7.2|) becomes 

Pn{S) = h{n) ■ det [es^+^-s,)'l^^^Q ■ 

Now use the argument of Theorem \A.2\ together with the identification of the normal- 
izing constant h{n) = l/e(l)"' in Theorem 17. 5[ □ 

Remarks. 

1. Another approach to Theorem 17.51 is via the theory of conditional L-ensembles in 
Borodin-Rains |i. Prop. 1.2]. 



2. In the translation invariant case e{i,j) = e{j — i), Stanley 67|, p. 90, Ex. 14] shows 
that 

Pniin — 1]) = chance that all sites are occupied 

is the coefficient of in the power series h{n)e{—z)^^ . To prove this using 
symmetric function theory, note from Corollary 17.71 that 



Pn{[n-l]) = hin) 



Here the Schur function S(„) = hn (where hn is the nth complete homogeneous 
symmetric function and h{n) is the normalizing constant in (|7.2p ). If E{z) = 
Xlr^o ^rZ^ is the generating function for elementary symmetric functions and 
H{z) = YT=ohrz'', Macdonald HH, p. 21, (2.6)] shows that H{z)E{-z) = 1. 
This gives Stanley's formula since in Corollary 17. 7^ the value of is e(r). 
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Example 26. (Descents in a random sequence) From Fact 5 of lSection 21 we are 
in the Toeplitz case witli e{j) = ("'^^7^) , e(z) = {1 — z)^^. Applying Corollary 17. 61 yields 
Theorem 13.11 

Example 27. (Descents in a uniform permutation) From MacMahon's formula 
(jS.ip . we are again in the Toeplitz case with e(j) = l/j\,e{z) = e^. As in Example O 
one can replace z by —z without changing determinants or the correlations functions. 
This proves Theorem 15.11 

Example 28. (Descents in a non-uniform permutation) From Stanley's formula 
(part a of Proposition 15. 2p . we are again in the Toeplitz case with e{j) = 1/jlq {q- 
factorial) . An identity of Euler allows one to write 

^(^) = n l-z(l- q)qm 

m>0 ^ ^'^ 

when 0<g<l,|2:| <1, but the elementary description of the correlation functions 
given in Proposition 15.21 seems more useful. 
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